February 5, 2008 



Interferometric multi-wavelength (sub)millimeter continuum study 
of the young high-mass protocluster IRAS 05358+3543 

H. Beuther 1 , S. Leurini 2 - 3 , P. Schilke 2 , F. Wyrowski, K.M. Menten 2 and Q. Zhang 4 

1 Max-Planck-Institute for Astronomy, Konigstuhl 17, 691 17 Heidelberg, Germany 
e-mail: beuther@mpia.de 

2 Max-Planck-Institute for Radioastronomy, Auf dem Hiigel 69, 53121 Bonn, Germany 
e-mail: name@mpifr-bonn.mpg.de 

European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany 
e-mail: sleurini@eso.org 
4 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
e-mail: qzhang@cfa.harvard.edu 



ABSTRACT 

Aims. We study the small-scale structure of massive star-forming regions through interferometric observations in several (sub)mm 
wavelength bands. These observations resolve multiple sources, yield mass and column density estimates, and give information about 
the density profiles as well as the dust and temperature properties. 

Methods. We observed the young massive star-forming region IRAS 05358+3543 at high-spatial resolution in the continuum emission 
at 3.1 and 1.2 mm with the Plateau de Bure Interferometer, and at 875 and 438 yum with the Submillimeter Array. The observations 
are accompanied by VLA 3.6 cm archival continuum data. 

Results. We resolve at least four continuum sub-sources that are likely of protostellar nature. Two of them are potentially part of 
a proto-binary system with a projected separation of 1700 AU. Additional (sub)mm continuum peaks are not necessarily harboring 
protostars but maybe caused by the multiple molecular outflows. The spectral energy distributions (SEDs) of the sub-sources show 
several features. The main power house mml, which is associated with CH 3 OH maser emission, a hypercompact Hn region and 
a mid-infrared source, exhibits a typical SED with a free-free emission component at cm and long mm wavelengths and a cold 
dust component in the (sub)mm part of the spectrum (spectral index between 1.2mm and 438yum a ~ 3.6). The free-free emission 
corresponds to a Lyman continuum flux of an embedded 13 M Q B 1 star. The coldest source of the region, mm3, has a ~ 3.7 between 
1.2 mm and 875 /jm, but has lower than expected fluxes in the shorter wavelength 438 yum band. This turnover of the Planck-function 
sets an upper limit on the dust temperature of mm3 of approximately 20 K. The uv-data analysis of the density structure of individual 
sub-cores reveals distributions with power-law indices between 1.5 and 2. This resembles the density distributions of the larger-scale 
cluster-forming clump as well as those from typical low-mass cores. The files for figures 2, 3, 4 & 6 are also available in electronic 
form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+AJ. 
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1. Introduction 

Single-dish (sub)mm continuum imaging in high-mass star- 
forming regions revealed various interesting properties of the 
dust and gas clumps harboring the star-forming clusters. The 
derived physical parameters include the gas masses, the gas 
column densities, the spectral indices and resulting dust prop- 
erties as well as the density dist ri butions of the gas on 
cluster scales flHaJcJiell et all l2000t Ivan der Tak et al l2000t 
Molinari et al.ll2000t iMueller et al.ll2002t iBeuther et al.l l2002bk 
Hatchell & van der Takll2003t I Walsh et alj|2003|: I Williams et all 
2004; Fontani et al. 2005; Williams etaTTl2005; Hill et al. 2005; 



Klein et all l2005t iFaundez et all l2004t iBeltran et all 120061) . 
However, all these studies had insufficient spatial resolution (at 
best 16000 AU at 2kpc distance) to resolve sub-structure within 
the star-forming regions, and thus were inadequate to study the 
properties of individual cluster members or unresolved multi- 
ple systems. A further limitation was that many of such stud- 
ies only covered one frequency band (870yum or 1.2 mm) and 



hence could not investigate the spectral index of the gas emis- 
sion, which allows to set constraints on the dust properties. 

To overcome some of these limitations we conducted a com- 
prehensive program observing the very young massive star- 
forming region IRAS 05358+3543 at high spatial resolution with 
the Submillimeter Array (SMAQ) in the 438 and 875 yt/m band, 
and with the IRAM Plateau de Bure Interferometer (PdBfl) in 
the 1.2 and 3.1 mm bands. These observations resolve the high- 
mass star-forming cluster into various sub-components at four 
different wavelengths, allowing to study the physical properties 
of individual cluster members or still unresolved, likely gravita- 
tionally bound multiple systems in great detail. 

IRAS 05358+3543 is part of a sample of 69 High-Mass 
Protostellar Objects (HMPOs) investigated in depth over the last 
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few years bvlSridharan eTail (120021): [Beuther et alJ (l2002bldldh : 
Willi ams et al.l d2004 120051) : lFuller et al.l d2005l) . At a kinematic 
distance of 1.8 kpc, the bolometric luminosity of the region is 
10 3 ' 8 L Q , and the integrated gas mass is ~300M(tJ3- H2O and 
Class 11 C H3OH maser emission was studied at high spatial 
resolut ion dTofaniet al.lll995t; iMinier et al.ll2000t IBeuther et all 
2002d|), but so far no cm continuum emission had been de- 
tected down to the threshold of 1 mJy, indicative of the the 
very young evolutionary stage of the region. Furthermore, high- 
spatial-resolution interferometric observations of the region re- 
vealed at least three molecular outflows, one of them being 
the most collimated massive molecular outflow known today 
(IBeuther et al.l l2002al) . These observations also resolved the 
evolving protocluster at 2.6 mm wavelength with a 4" x 3" 
beam into three sub-sources (mml to mm3) with gas masses 
between 37 and 50M Q 3 . Fig. Q] gives an overview of the large 
and sm all-scale dust continuu m emission from previous obser- 
vations (Beuther et al. 2002b a]). Two of the three identified out- 
flows emanate from the vicinity of mml. McCaughrean et al. (in 
prep.) detected with Keck observations two mid-infrared sources 
in the region, one being spatially coincident with the main mm 
continuum peak mml. This mid-i nfrared source was recently 
resolved by Long more et al.l d2006l) into a double-source sep- 
arated by 1" corresponding to 1800AU. Near-infrared polari- 
metric imaging also revealed a d eeply embedded source associ- 
ated with mml dYao et al.ll2.000l) . Previously identified infrared 
sources, some of them were suggested to be the outflow driving 
sources (Porra s et al. 2000), were shown to be offset from the 
main mm and mid-infrared source s and hence unlikely to be the 
main power houses in the region (Beut her et alj |2002a). 

In this paper, we present the (sub)mm continuum observa- 
tions from all four wavelengths bands. The corresponding spec- 
tral line data will be published separately by Leurini et al. (in 
prep.). For the nomenclature of this paper, we use the term proto- 
star (and corresponding proto-cluster, proto-binary) for massive 
young stellar objects that are still actively accreting, independent 
of whether they have started hydrogen burning already or not. 



2. Observations and Data Reduction 

We observed IRAS 05358+3543 with three interferometers 
(SMA, PdBI and VLA) in ten different configurations or spectral 
setups. The most important observational parameters are given 
in TableQ] The velocity of rest was vi sr = — 17.6kms _1 , and the 
flux accuracy of the observations is estimated to be good within 
20%. Below we discuss a few additional details to the different 
observations. 



2. 1. Submillimeter Array Observations at 875 //m 

The SMA observations in the compact and extended configura- 
tion were centered at 348/338 GHz (upper and lower sideband, 
USB/LSB, 875 yum). The projected baselines ranged between 15 
and 252 H. The short baseline cutoff implies that source struc- 
tures > 16" are filtered out by the observations. The phase 
center of the observations was R.A.(J2000.0) 05 A 39 m 13 s .07 
and decl.(J2000.0) +35°45'51".2. The zenith opacities, mea- 
sured with the NRAO tipping radiometer located at the Caltech 
Submillimeter Observatory, were good during both tracks with 
t(348GHz) ~ 0.18 (scaled from the 225 GHz measurement via 



3 The cited gas masses are at a 50% level of the masses originally 
presented in IBeuther et al.l ( 2002b a) because of different assumptions 
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Fig. 1. Overview of the region. The large emission shown in 
grey-scale with dotted contours is the 1 .2 mm continuum emis- 
sion o bserved with the MA MBO array at the IRAM 30m tele- 
scope ( Beuther et al. 2002b). The full contours (enlarged in the 
inlay-box) show the 2.6 mm continuum data from the PdBI 
at a spatial resolution of 4.0" x 3 .0" with the labeled three 
identified mm continuum sources (Beut her et ail l2002al) . The 
1.2 mm MAMBO data are contoured from 106 mJy beam" 1 
to 954mJybeairT 1 (step 106 mJy beam" 1 ), the 2.6 mm PdBI 
are contoured from 6.6mJybeam" 1 to 22.0mJy beam" 1 (step 
2.2 mJy beam" 1 ). The (0/0) position of the large map is the posi- 
tion of the IRAS 05358+3543 source which is derived from the 
IRAS 12yum image, whereas the mm continuum peak correlates 
with the IRAS 100 yum peak. The (0/0) po sition of the in l ay is th e 
phase center of the older observations bv lBeuther et al.l (|2()02a). 



r(348GHz) ~ 2.8 x r(225GHz)). The correlator had a bandwidth 
of 1.968 GHz and the channel separation was 0.8 125 MHz. 
Measured double-sideband system temperatures corrected to the 
top of the atmosphere were between 150 and 500 K, mainly de- 
pending on the elevation of the source. The initial flagging and 
calibration was done with the IDL superset MIR originally de- 
velop ed for the Owens Valley Radio Observatory ( Scovill e et al.l 
1 19931) and adapted for the SMA0 The imagin g and data analysis 
was conducted in MIRIAD dSault et aljl995l) . During the obser- 
vations on November 1 1th, 2004, the position of the primary cal- 
ibrator 3C1 1 1 was wrong in the catalog by ~0.6". Therefore, we 
self-calibrated our secondary phase calibrator 0552+398 (16.3° 
from the source) shifting it in the map to the correct position. The 
solutions were then applied to IRAS 05358+3543. The 875 jum 
continuum data were produced by averaging the apparently line- 
free part of the USB and LSB data. The corresponding line emis- 



for the dust opacities ( IBeuther et all2005h 



4 The MIR cookbook by Charlie Qi can be found at 
http://cfa-www.harvard.edu/~cqi/mircook.html 



sion will be discussed in a separate paper by Leurini et al. (in 
prep.). 

2.2. Submillimeter Array Observations at 438 /jm 

The weather conditions for the 438/vm/690GHz observations 
were excellent with a zenith opacity t(230GHz) between 0.03 
and 0.05 throughout the night. This corresponds to zenith opac- 
ities at 690 GHz be tween 6 and 1.0 (t(690GHz) ~ 20 x 
(t(230GHz) - 0.01), lMassonlll994l) . The phase center and cor- 
relator setup were the same as for the 875 //m observations. The 
primary beam of the SMA at these high frequencies is ~ 18". 
The covered frequency ranges were 679.78 to 681.75 GHz and 
689.78 to 691.75 GHz. Measured double-sideband system tem- 
peratures corrected to the top of the atmosphere were between 
1350 and 10000K, depending on the elevation of the source and 
the opacity. The sensitivity was limited by the strong side-lobes 
of the strongest emission peaks. Since the third main mm con- 
tinuum source is approximately at the edge of the SMA primary 
beam at this frequency, we applied a correction for the primary 
beam attenuation to the final image to determine the fluxes more 
properly. Although the synthesized beam of the continuum im- 
age is 1.4" x 0.9", we work only on a smoothed version with a 
1 .9" x 1 .4" beam, comparable to the 3 mm data. 

2.3. Plateau de Bure Interferometer Observations at 1.2 and 
3.1 mm 

The PdBI was used in two different frequency setups in 2003 
and 2005. A first frequency setup was performed in four tracks 
between January and October 2003 in the B-C and D con- 
figurations of the array. The reference center is R.A. (J2000) 
05*39 m 13 s .O and decl. (J2000) 35°45'54". The 3 mm receivers 
were used in single sideband mode and tuned to 96.6 GHz; 
the 1 mm receivers, in double sideband mode, were tuned to 
241.85 GHz, (USB). One 320 MHz unit was placed to obtain a 
continuum measurement at 3.1 mm. Four overlapping 160 MHz 
units were used to obtain the 1.2 mm continuum measurement. 
The system temperatures varied between 100 and 150 K in the 

3.1 mm band, and between 300 and 600 K in the 1.2 mm band. 
The second frequency setup was observed in 2005 in the A and 
B configurations. The 1mm receivers were tuned in lower side- 
band to 241.2 GHz. Three 160 MHz units were used to get the 
continuum measurement. Measured system temperatures in the 

1.2 mm receivers ranged between 240 and 1000K, with one re- 
ceiver measuring system temperatures of 900 K. The data cali- 
bration were performed in CLIC, and the imaging in MAPPING, 
both packages are part of the GILDAS software suit^B 

2.4. Very Large Array archival data at 3.6 cm 

Searching the Archive of the Very Large Array (VLA), we found 
A+B array data observed at 3.6 cm wavelength in 2002 and 2003 
(project AR482). The data were reduced with the standard AIPS 
procedures. 

3. Results 

3. 1 . Small-scale dust/gas distribution 

The high-spatial-resolution multi-frequency data allow to dis- 
entangle considerable sub-structure in this massive young star- 

5 See http:/ /www.iram.fr/IRAMFR/GILl5ASl 



forming cluster. Fig. [2] presents overlays of the 1.2 mm and the 
875 //m dust continuum emission. The top-panel shows images 
at both wavelengths with comparable resolution (BCD config- 
uration with the PdBI at 1 .2 mm, and extended+compact con- 
figuration with the SMA at 875 jum). The 3.1mm and 438 jum 
images with slightly lower spatial resolution will be discussed 
below ( ^3.2t . We clearly resolve many sub-features compared to 
the more simple three source s previously detecte d at lower spa- 
tial resolution (see Fig.fTland lBeuther et al.l20 02a). Furthermore, 
we have additional 1 .2 mm continuum data with even higher spa- 
tial resolution (0.6" x 0.4", AB configuration with the PdBI) 
shown in the bottom-panel of Fig. [2] While this highest-spatial- 
resolution dataset filters out nearly all large-scale emission it 
clearly identifies the most compact sources. We could have com- 
bined all four configurations, but since the most extended con- 
figurations dominate this image, we refrained from that. 

The main (sub)mm continuum source mml, which i s at the 
center of two molecular outflows dBeuther et al.ll200"2al) . is re- 
solved into two separate (sub)mm continuum peaks, mm la and 
mmlb. This is already observable at ~1" in Fig. [2] (top panel) 
but more clearly discernable in the highest-resolution data in the 
bottom panel of Fig. [2] and the zoom in Fig. [3] The compact 
3.6 cm continuum source and the Class n CH3OH maser are as- 
sociated with mml a. The cm continuum source is unresolved at 
the given spatial resolution with a pea k flux of 0.62mJybe am~ 1 
(below the detection limit of 1 mJy by Sridh aran et al.l2002h . and 
we only show the 50% contou r level in Figs.[2]&|3] The closest 
H2O maser feature detected bv lTofani et al J d 1 995b is — 1" north- 
east of mm la and not clearly associated with one or the other 
peak as well. It should be note d that this feature was not de- 
tected by Beut her et al.1 d2002dl) . probably due to variability of 
the feature. The projected separation between mmla and mmlb 
amounts to ~0.94", corresponding to a projected linear sepa- 
ration of ~1700AU at the given distance (1.8kp c ). The main 
mid-infrared source detected by Longmor e et al.l (|2()06) coin- 
cides well with the position of mmla (Fig. [3}. However, while 
the secondary source mmlb is ~1700 AU toward the south-east, 
the secondary mid-infrared source is ~1800AU south-west of 
mmla. Hence, mmlb and the secondary mid-infrared source are 
apparently independent indicating even more sub-structure in the 
evolving young cluster. 

The previously identified source mm2 is resolved even more, 
at least into two sub-sources mm2a and mm2b (Fig.[2j. However, 
while mm2a shows the same peak position in the 1 .2 mm band 
and the 875 /im band, the peaks associated with mm2b are sepa- 
rated by approximately 0.4". The mm peak mm2c is identified in 
the 1.2 mm image but not at shorter wavelengths. Inspecting the 
longer wavelength 3.1 mm data in Figure [4] one finds an elon- 
gation in this lower-spatial-resolution dataset toward the mm2c 
position as well (labeled "jet?" in Fig. 0). Since also the spec- 
tral index between the 1 .2 mm and 875 //m toward this position 
is very low (of the order 1.5, see 93.4l and corresponding Fig. |6j 
this feature may be due to an underlying outflow/jet. The fea- 
ture mm2d in Fig. [2] is ag ain discernable in both bands. One of 
the H2O maser features bv lTofani et al.l d 1995b is associated with 
this latter (sub)mm continuum emission feature. 

The third source mm3 remains a single compact object even 
at the highest spatial resolution. The lower-resolution 3.1mm 
image in Fig. [T] (inlay) shows an elongation from mm3 toward 
the south-west, and the higher-resolution images in Fig.|2]show 
weak emission there as well. Since these weak emission fea- 
tures show different peak positions in the 1 .2 mm and the 875 /im 
bands, we refrain from further discussion of this feature. 



Table 1. Main observation parameters 
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The Table shows the observing frequency A, the observatories, the dates of observation, the array configurations, the calibrators, the resulting 
synthesized beam and the lcr value of the final maps. 



Protostellar condensations should be detectable in all 
(sub)mm wavelengths bands and not be filtered out in very ex- 
tended configurations. Therefore, we identify at ~0.5" spatial 
resolution (~ 900 AU) four compact protostellar sources in the 
region (mmla, mmlb, mm2a, mm3) , the positions are listed in 
Table El 

Table 2. Positions of the four identified protostellar sources 



Source 


R.A. [J2000] 


Dec. [J2000] 


mmla 


05:39:13.08 


35:45:51.3 


mmlb 


05:39:13.13 


35:45:50.8 


mm2a 


05:39:12.76 


35:45:51.3 


mm3 


05:39:12.50 


35:45:54.9 



3.2. The multi wavelength dataset 

Since our 3.1mm data are of lower spatial resolution, we can 
smooth the higher frequency data to this resolution (1.9" x 
1.4") and investigate the SEDs of the relatively larger sub- 
components. Figure H] presents the compilation of all four con- 
tinuum images at the spatial resolution of the 3.1 mm continuum 
dataset. We clearly identify three (sub)mm continuum peaks at 
all wavelengths, the formerly known sources mml to mm3. The 
vertical lines in Figure @] locate the R.A. positions in the three 
lower-frequency bands (3.1 mm, 1 .2 mm and 875 /mi). While the 
positions of mml and mm3 coincide very well between the dif- 
ferent bands, the position of mm2 shifts progressively toward 
the south-east going from 3.1 mm to 875 /im. However, compar- 
ing the 1 .2 mm and 875 fim lower-resolution images in Figure 
@] with the higher-spatial-resolution images presented in Figure 
|2] we see that the source mm2 consists of at least two sub- 
sources, mm2a and mm2b. While the 1 .2 mm and 875 /mi po- 
sitions of mm2a coincide well (Fig.|4]and full circles in Fig. |4j, 
the peak position of mm2b is shifted toward the south west 
going from 1.2 mm to the 875 /mi data. Based on the lower- 
resolution 3.1 mm image, it appears likely that at higher spatial 



resolution mm2b would be even closer to mm2a going to even 
longer wavelengths. This positional shift of mm2b as well as 
the non-detection of mm2c in the 875 /im indicates that the dust 
emission at these positions is not of typical protostellar nature, 
but caused by other processes, e.g., the protostellar outflows in 
the region. As already mentioned in 33. 11 the 3.1 mm image ex- 
hibits a smaller peak between mml and mm2 which is approxi- 
mately associated with mm2c in the shorter-wavelength higher- 
resolution data (Fig.f2]l- Taking additionally into account the low 
spectral index toward this position (see belo w), this feature i s 
likely produced by an underlying jet/outflow dRevnoldsl fl986). 
hence we labeled it the jet-feature in the 3.1 mm image. 

The 438 /im image requires some additional explanations 
because interferometric imaging at these frequency is still ex- 
ploratory and very difficult. The sub-source mml is at the 
phase center and the most reliable feature. The emission peak 
of mm2 does not shift further to the south-east as maybe ex- 
pected from the lower-frequency data, but the R.A. position co- 
incides approximately with the 1.2 mm peak position whereas 
the Declination position is even closer to the 3.1 mm position. 
Furthermore, the emission around mm2 in this band extends a 
little bit further to the north than in the other images. We are 
not entirely certain whether this is an artifact of the compara- 
bly lower data quality in this band, or whether it may even be 
real because many of the molecular line maps (e.g., CH3OH, 
H2CS or C 34 S) have an additional emission peak north of mm2 
as well (Leurini et al. in prep.). Since high-spatial-resolution im- 
ages at these wavelengths are extremely ra re so far (only 2 are 
published to our knowledge, Or ion-KL by Beut her et alj |2006b 
and TW Hydra bv lOi et al.ll2006l) . it is hard to set additional con- 
straints here. The location of mm3 is right at the edge of the 
primary beam of our dataset (the phase center was chosen like 
that because we were also interested in the CO(6-5) emission 
from the main outflow emanating from mml), and the measured 
fluxes toward mm3 are barely above the noise (~2.5<x). To at 
least derive an upper limit for the emission from mm3, we cor- 
rected the image shown in Figure [4] for the primary beam atten- 
uation. The derived fluxes shown in Tableland used below for 
the mass/column density estimates, and spectral energy distribu- 
tions are measured from the primary-beam-corrected image. 
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Fig. 2. (Sub)mm continuum images toward IRAS 05358+3543 at 1.2 mm (PdBI) and at 875//m (SMA). The top panel shows in 
grey-scale with dotted contours the 1.2 mm data (BCD configuration) and in full contours with increasing thickness the 875 jum 
data. The contouring is done from 10 to 90% of the peak emission of each image (step 10%), with peak flux values of 58.9 and 
164.3 mJy beam -1 for the 1 .2 mm and 875 pm bands, respectively. The bottom-panel shows in grey-scale with dotted contours again 
the same 1 .2 mm data (BCD configuration), and in full contours the 1 .2 mm data observed at even higher spatial resolution (0.6" X 
0.4") in the AB configuration. The contouring of the latter dataset is from 20 to 90% (step 10%) of the peak flux 29.8 mJy beam -1 . In 
both panels, the white ellipse s hows the VLA 3.6 c m 50% contour level (0.3 1 mJy). The st ar, triangles and sq uares mark the positions 
of th e Class n CH3OH ma ser dMinier et all 2000). the H2O masers (open triangles from lTofani et al.l ll995 and small filled triangle 
from lBeuther et al1l2002dl) and the mid-infrared sources (McCaughrean et al. in prep.). The newly identified double-source mmla 
and mmlb is labeled as well. The box outlines the region shown in Figure[3] The files for this figure are also available in electronic 
form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr ( 130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat? 'J/A+A/. 



3.3. Mass and column density estimates 

Assuming that the (sub)mm continuum emission is optically 
thin and thus tracing the whole dust column densities with 
a gas to dust mass ratio of 100, we can derive t he gas col- 
umn d ensit ies Afa, and gas mas ses M following iHildebrandl 
d 1983b and iBeuther et alj d2002bl) . The following sections will 
discuss the spectral indices toward the various sub-regions in 



IRAS 05358+3543 in more detail, but for the mass calculations 
we adopt a dust opacity index ft of 1.5, corresponding to dust 
opacities k of 0.2, 0.9, 1.4, and 4.0cmV at 3100, 1200, 875, 
and 438 pm, respectively. Although the region contains gas com- 
ponents at various temperatures (see Leurini et al. in prep.), we 
assume an average temperature of 50 K for the mass estimates. 
Only for mm3, we use a lower temperature of 20 K as derived 
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Fig. 3. Zoom into the central region of mmla and mmlb as 
shown in Figure [2] The grey-scale with dotted contours shows 
again the same 1.2 mm data (BCD configuration), and in full 
contours the 1.2 mm data observed at higher spatial resolution 
(0.6"x0.4") in the AB configuration. The contouring is the same 
as in Figure [2] The white ellipse shows the VLA 3.6 cm 50% 
contour level (0.31 mJy). The star, triangle s and squares mar k 
the positions of the Class n CH3OH maser dMinier et al.l l2000). 
the H2O masers dTofani et al.l l995) and the mid-infrared sources 
(McCaughrean et al. in prep.). The black-and-w hite circles show 
the do uble mid-infrared sources discussed by Longmore et al.l 
(2006!). The 0.5" positional difference between the main mid- 
infrared peak positions by McCaughrean et al. and Longmore 
et al. can be accounted for by relatively poor pointing accura- 
cies and the difficulties of registering mid-infrared data. The files 
for this figure are also available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+A/ 



from the spectral energy distribution shown in 33.41 Given the 
uncertainties of the dust opacities and the temperatures, we es- 
timate the masses and column densities to be accurate within a 
factor 4. Table[3]presents the derived masses for the various sub- 
sources estimated from the multi-wavelength dataset. 

As expected, the derived masses are all below the previous 
single-dish 1.2 mm and lower-resolution 3 mm PdBI mass esti- 
mates (©. Some part of the discrepancy can be solved by the 
different assumptio n for the dust opacity J 3 used in the previ- 
ous studies (J3 = 2. iBeuther et alj|2~002blal) and here (J3 = 1.5). 
However, the difference in /3 accounts to a mass difference of 
approximately a factor 2 at 1.2 mm and 875 pm, and the larger 
effect is the missing flux in the new interferometer data. The in- 
tegrated single-dish fluxes at 1 .2 mm and 850 pm are 6 Jy and 
28.7 Jy (Willi ams et alj [2004b . respectively, whereas the inte- 
grated fluxes in our high-spatial-resolution 1.2 mm (BCD con- 
figuration) and 875 pm images (Fig. [2]) are 0.85 Jy and 1.37Jy, 
respectively. Thus, we are recovering only between 15% and 5% 
of the flux - and hence mass - in these interferometric images. 
The difference in flux recovery can already be recognized by 
just comparing the values for the 1.2 mm data in the AB and 
BCD configurations, respectively (Table (3). Therefore, we omit 
the mass estimates of the higher-resolution data in Table [3] 



Table 3. (Sub)mm continuum fluxes and derived masses/column 
densities 



Source 5i„ t . 


S Peak 


M' AT* 


[mJy] 


r mJy , 
L beam J 


[M s ] [cm- 2 ] 


3.1 mm (AB) with 1.9" x 1.4" beam 


mml 20 


10.3 


10.6 1.7 10 i4 


mm2 10 


5.9 


5.3 9.7 1 23 


mm3 d 13 


6.5 


18.5 2.9 1 24 


1.2 mm (BCD) with 1.9" x 1.4" beam 


mml 215 


91.8 


5.0 6.6 10 23 


mm2 64 


38.3 


1.5 2.7 1 23 


mm3<' 39 


37.6 


2.7 8.1 10 23 


875 pm with 1.9" x 1.4" beam 


mml 529 


270.8 


3.8 6.0 1 23 


mm2 199 


110.1 


1.4 2.4 10 23 


mm3 rf 165 


137.5 


3.9 1.0 10 24 


438 pm with 1.9" x 1.4" beam" 


mml 8630 


3589 


6.5 8.4 1 23 


mm2 1 1450 


3749 


8.6 8.8 1 23 


mmS^ 941 


941 


3.2 9.9 1 23 


1.2 mm (AB) with 0.60" x 0.44" beam 


mmla 42 


29.8 


1.0 2.2 1 24 


mmlb 26 


22.3 


0.6 1.6 10 24 


mm2a 3 1 


18.3 


0.7 1.3 10 24 


mm3 rf 27 


22.8 


1.9 4.9 1 24 


1.2 mm (BCD) with 1.26" 


x 0.84" beam 


mmla 67 


58.9 


1.1 10 24 


mmlb 62 


54.8 


9.8 10 23 


mm2a 67 


39.9 


7.2 10 23 


mm2b 33 


33.3 


6.0 10 23 


mm2c 65 


32.1 


5.8 10 23 


mm2d 41 


26.5 


4.8 10 23 


mm3 d 103 


48.0 


2.6 10 24 


875 pm with 1.14" x 0.57" beam 


mmla 225 


164.3 


1.5 10 24 


mmlb 124 


106.4 


9.6 10 23 


mm2a 97 


62.5 


5.7 10 23 


mm2b 69 


58.6 


5.3 10 23 


mm2c V 






mm2d 68 


46.2 


4.2 10 23 


mm3<' 212 


116.1 


3.5 10 24 



" No clear emission peak in that band. 

* Fluxes in this band are measured from the primary-beam corrected 
image. 

c The flux values are only 2.5<r. 

d A colder temperature of 20 K was used for mm3 ( 33. 4t , 

' The masses and column densities are estimated to be accurate within 

a factor 4. 



An additional effect arises from the various wavelength ob- 
servations: Comparing the three mass estimates derived from 
3.1 mm, 1.2 mm and 875 pm data with the same spatial resolu- 
tion (1.9" x 1.4", Table O, the masses at 1.2 mm and 875 pm 
are approximately the same, whereas the masses obtained from 
the 3.1 mm data are more than a factor 2 higher. We will discuss 
the spectral energy distributions in detail below, but already this 
comparison indicates that at 3 mm additional non-dust contribu- 
tions to the continuum flux have to be taken into account. 

Regarding the flux measurements and mass/column density 
estimates in the 438 pm band, the values for mml are reason- 
able compared to the other bands. Since the 438 pm flux for 
mm3 is only 2.5<x limit, the derived mass and column density 
can be considered as an upper limit. The strangest values are 
for mm2 in the 438 pm band because they are higher than for 
mml whereas this is always the other way round in the longer 



wavelength bands. As discussed in the previous section, with the 
current data is is difficult to decide whether this is an artifact 
of the dataset, or whether there are additional contributions to 
the continuum emission in this band as maybe suggested by the 
strong line emission peaks a small distance to the north of mm2 
(Leurini et al. in prep.). 

3.4. Spectral energy analysis from the multi-wavelength 
(sub)mm continuum images 

In spite of the complexity of the region, we can derive the SEDs 
from the lower-resolution images toward mml, mm2 and mm3 
(Fig.H]). The emission in the 3.1 mm and 438 yum bands appears a 
little more extended which may be due to a better sampling of the 
uv-plane at the corresponding scales. Therefore, we refrain from 
using integrated fluxes for the various positions but rather use the 
peak fluxes listed in Table [3] The peak fluxes are derived at the 
peak position of each image. Figure[5]presents the final (sub)mm 
SEDs for the three positions at the given spatial resolution. 

For all three (sub)mm peaks mml to mm3, the 3.1 mm point 
is consistently too high to be fitted with a single power law which 
is expected for pure dust continuum emission at (sub)mm wave- 
length. Therefore, the (sub)mm power-law fits at all three posi- 
tions were derived from the shorter wavelength data only. 

Source mml: The 1.2mm to 438 yum data are well fitted by 
a single power-law S oc y" with a ~ 3.6 + 0.3. This power- 
law exponent corresponds in the Rayleigh-Jeans regime to a 
dust opacity index j3 = a — 2 of 1.6 ± 0.3, observed often in 
high- and low-mass star-forming regions (e.g. lAndre et alj2000t 
Beut her et al.ll2006ah . The low-frequency 8.3 GHz data from the 
VLA are likely from a hypercompact Hn region. Therefore, we 
can fit the whole spectrum from 3.6 cm to 438//m with a uni- 
form density hypercompact Hn region at low to intermediate fre- 
quencies and a dust component in the short mm to submm range 
(Fig. [5] top-left panel). This is a typical spectrum one expects 
from a young massive protostar where the c entral source has 
alread y started hydrogen burning. Following Rohlf s~& Wilson! 
(2006), we can calculate the emission measure (EM) using the 
approximate turnover frequency of 35 GHz and en electron tem- 
perature of ~ 10 4 K. The derived EM of ~ 5.4 x 10 9 cirT 6 pc is 
very high, about one to two orders of magnitude high er than val- 
ues fo und in typical ultracompact Hn regions (e.g. jKurtz et alj 
1 1 9941) . Taking the optical thin free-free fl u x of ~6.4mJy at 
60 GHz (Fig. O and following iKurtz et alj d 19941) . we calcu- 
late a Lyman continuum flux of ~ 10 45 4 photons s , which cor- 
responds approximately to a B 1 Zero- Age-Main-Sequence star 
with a luminosity of ~ 10 3 72 L Q and a stellar mass of ~ 13 M 
dPanagial 19731 : ILangl 19921) . Based on the strong, collimated out- 
flow that emanates from the vicinity of mml with and outflow 
rate of 6 x 10 4 M Q an d a dynamical timescale of 3.6 x 10 4 yrs 
dBeufher et alj|200"2al) . it is likely that this source is still in its 
accretion process, and therefore, we still call it a protostar. The 
luminosity derived from the free-free emission is of the same or- 
der a s the total luminosity (~ 10 3 8 L o ) derived from the IRAS 
data ( Srid haran et alj|2002l) . indicating that the luminosity of the 
whole region is really dominated by the most massive protostar. 

Source mm3: We can fit the 1 .2 mm and 875 /mi data reason- 
ably well with a power-law with a ~ 3.7 + 0.3 (corresponding to 
/? ~ 1 .7 + 0.3). The 3.1 mm fluxes are too high for that fit. Taking 
our 8.3 GHz VLA 3<x upper limits of 56 fiJy, it is again possible 
to fit the low-frequency part of the spectrum with a hypercom- 
pact Hn (Fig.|5]bottom-left panel). A point of caution here is that 
the turnover frequency from the optically thick to the optically 
thin regime is relatively high, between 70 and 100 GHz. One 



possible way to explain this is that the region is so young that the 
forming HCHu region is still quenched and thus optically thick 
up to relatively high frequencies without having the time and 
power to heat the surrounding gas. Of particular interest is the 
high-frequency data-point at 438 fim which is below the expec- 
tation of the power-law fit from the 1.2 mm and 875 fim. Since 
we have no cm nor compact line detection toward mm3 (see 
Beutheretal]|2002 a|and Leurini et al. in prep.), this sub-source 
is a potential candidate for a still very cold massive core at even 
earlier evolutionary stages. For such a cold core, the peak of the 
SED moves to longer wavelength and one does not expect a nor- 
mal Rayleigh- Jeans-like power-law at short submm wavelengths 
anymore. Therefore we tried to fit the high-frequency SED with 
grey-body functions with varying temperature (blackbody func- 
tions modified by 1 - e Tv taking into account the finite optical 
depth t y of the dust at all frequencies). The two bottom-panels 
in Fig. [5] show such grey-body functions with temperatures be- 
tween 50 and 15 K (we also varied [} between 1.7 and 2). One 
has to keep in mind that mm3 is at the edge of the primary beam, 
that the imaging is difficult and that the data are primary-beam 
corrected. Nevertheless, these grey-body functions indicate that 
most of the envelope dust and gas is likely at < 20 K. 

Source mm2: As already mentioned in the previous sections, 
mm2 is the strangest of the three sources. We do not find any rea- 
sonable fit for the 1.2 mm to 438 yum. The fits presented in Fig. [5] 
(top-right panel) show the extreme power-laws fitting the 1 .2 mm 
and the 438 fim data. The 875 yum data point is below these fits. 
As pointed out previously, it may be possible that the unpropor- 
tionally high 438 /im flux density could be partly an artifact of 
this dataset. However, there may also be additional contributions 
in this band (e.g., from the adjacent molecular line peak) which 
have to be taken into account separately. Furthermore, the data 
used for these SEDs are of lower spatial resolution smearing out 
mm2a and mm2b. Since mm2b may be caused by an outflow/jet 
with a lower spectral index (see the following paragraphs, this 
will likely increase the complication of interpreting mm2's SED 
(see also the cases of LI 1 57 or IRAS 20293+3952, lGueth et all 
120031: iBeuther et al.ll2004l) . In principal, a flattening of the SED 
in the long wavelength regime could also be produced by a 
change of the dust grain size distribution towards larger fractions 
of large grains within a protostellar disk. However, this appears 
an unlikely scenario here because a low spectral i ndex is also ob- 
serve d in the larger-scale envelope of the region (Williams et al. 
2004). Based on these data and the complicated structure of the 
region, we are not able to draw further conclusions from the SED 
of mm2. 

To derive additional small-scale information about the 
(sub)mm dust continuum emission, we used the higher-spatial- 
resolution 1.2 mm and 875 yum data to produce a spectral index 
map of the (sub)mm continuum emission. Therefore, we restored 
both datasets with the same synthesized beam of 1.2" x 0.8". 
To account for absolute positional differences due to non-perfect 
phase calibration, we shifted the 875 /mi image by 0.18" to the 
east aligning the peak-positions of mm la and mm3. Then the 
ratio map was constructed from these two images (Figure |6). 
The spatial distributions of the two restored maps shown in the 
top-panel of Figure|6]show that the spatial scales traced in both 
bands are approximately the same, and thus a comparison of 
both images appears reasonable. Taking into account the flux un- 
certainties of 20%, the different uv-sampling and the positional 
shift, we estimate the uncertainties of this spectral index map of 
the order +0.5. Edge effects of rising spectral index toward the 
north-east of mmla or toward the south-west of mmlb are most 
likely an artifact. Similarly, the contrast of a low spectral index 




Fig. 5. Spectral energy distributions (SEDs) for mml, mm2 and mm3. For mml and mm3, the data are fitted with a uniform density 
Hn region at low frequencies and a power-law with exponent 3.6 and 3.7 in the (sub)mm part of the spectrum, respectively. For 
mm3, we present also various grey-body functions varying the temperature T and the dust emissivity index /? as labeled in the plot. 
The arrows in the mm3 SEDs mark the upper limits (lcr noise at 8.3 GHz, and the tentative 941 mJy detection corresponding to 
2.5cr at 685 GHz). For mm2, we show two extreme power-laws fitting the 1.2 mm and the 438 /-im but missing the 865 fim point. We 
measured the fluxes toward the peak positions of each image in Fig. [4] respectively. 



north of mm3 and a high spectral index south of mm3 are po- 
tential artifacts caused by the division of the different datasets. 
However, the variations of the spectral indices throughout the 
region which are not edge effects are likely to be real. 

The spectral indices we derive from Figure [6] toward the 
four protostellar condensations mmla, mmlb, mm2a, and mm3 
(Tabled are 3.6, 2.7, 2.3, and 3.6, respectively. The spectral in- 
dices toward mmla and mm3 - approximately the same as in the 
lower-resolution data discussed before - correspond t o a rather 
normal dust opacity index yS of 1.6 (e.g., Draine 2006), whereas 
the spectral indices toward mmlb and mm2a are significantly 
lower. Since we aligned the two wavelength images simultane- 
ously with the peak positions from mmla and mm3, the lower 
values toward the other sources are likely to be real and no ar- 
tifact from the shifting. The positional shifting of mm2b from 
1.2 mm to 875 /mi (Fig. |2|i is reflected in the ratio gradient to- 
ward that region. Since mm2c is not detected at 875 /mi at all, 
the corresponding spectral index is hence very low (of the order 



1.5) indicative of an underlying jet (lRevnoldslll986l) . Another 
interesting position appears to be west of mm2a, where particu- 
larly high values > 4 are found in the ratio map. Figure|6]clearly 
demonstrates the complexity of spectral index studies and the 
large spread in spectral indices over very small spatial regions. 
Possible physical implications will be discussed in §4.21 

3.5. Radial structure analysis from the multi-wavelength 
(sub)mm continuum uv-data 

To investigate the structure of the protostellar sub-sources de- 
tected in the continuum, the data were analyzed directly in the 
uv-domain avoiding the uv-sampling and cleaning difficulties. 
To decrease the influence of the other sub-sources on the analysis 
of a particular object, two-dimensional Gaussians were fit as a 
first approximation to the highest resolution images and the cor- 
responding models were removed from the uv-data, resulting in 
"clean" uv-datasets for each of the three main sub-sources. Then 



the data for each sub-source were binned in uv-annuli around 
the peak positions and averaged vectorially. The averaged am- 
plitudes versus uv-radius are shown for the 3 sources at the 3 
observed wavelengths in Fig. [7] Toward mml, we selected the 
the mm la peak position as the center, but obviously most base- 
lines trace structures combining mmla and mmlb, and only the 
largest baselines (>250k/l) are only due to mmla. Therefore, 
the derived intensity and density profile is for the multiple sys- 
tem mmla/mmlb. Most of the data can be fit reasonably well by 
straight lines in the log-log plots, showing that the radial struc- 
tures of the sources are close to power-law profiles. The data are 
inconsistent with Gaussians sources. Some of the data show ev- 
idence for constant amplitudes at the largest uv-radii, hinting at 
point source components at the smallest scales. Deviations from 
the straight lines are due to more complicated morphology im- 
printed on power-law profiles, e.g. the substructures in mml and 
mm2. 

Table [4] shows the results of the power-law fits to the data. 
Assuming that the emission is due to dust at all frequencies, the 
radial profiles should be the same at different frequencies. The 
differences that can be seen in the data therefore shows the ob- 
servational limitatio ns of this method of radial profile analysis. 
lLoonev et al. (2003) give a relation between the power-law pro- 
files in the image and uv-planes, as well as how the observed uv- 
power-law indices V(s) are composed by the density and temper- 
ature structure: V(s) ~ s p+q ^ , where p and q are the power-law 
indices of the density n ~ r~ p and temperature structure T ~ r~ q , 
respectively. With a typical temperature power-law index of 0.4, 
density indices from 1.5 to 2.0 translate to -1.1 to -0.6 in the uv 
plane. Inspecting the observed values in Table [4] shows that the 
millimeter cores have density profiles in this range, with more 
observed indices closer to the p — 1.5 than to the p = 2.0 pro- 
file. 

Table 4. Results of power-law fits to the continuum amplitudes 
versus uv-distances of the embedded protostellar sources 





mmla" 


mm2 




mm3 


a(3.1mm) 


-0.59 


-0.78 




-0.61 


a( 1.2mm) 


-0.86 


-0.71 




-0.53 


Q-(875/im) 


-0.93 


-1.05 




-0.94 


a av. 


-0.8 + 0.2 


-0.85 ±0.2 


-0.7 


±0.2 


p(3.1mm) 


2.0 


1.8 




2.0 


p{ 1.2mm) 


1.7 


1.9 




2.1 


p(875yum) 


1.7 


1.5 




1.7 


p av.* 


1.8 ±0.3 


1.75 ±0.3 


1.9 


±0.3 



" The fit was centered on mmla but it covers mainly structures combin- 
ing mmla and mmlb. 

* The error for the average p takes also into a account an error for q of 
±0.2. 



4. Discussion 

4.1. Protostellar Powerhouses and other sources 

The double-structure mmla and mmlb is separated by only 
~1700AU projected on the plane of the sky. This double- 
system may well be a massive bound binary system, similar 
to the the potential high-mass protocluster or b inary system 
W3(H 2 0) dWvrowski et : al.ll999tlChen et al.l2006l) . Considering 
also the third source d etected at mid-infrared wavelengths by 
lLongmore et al . ( 2006), there may be three sources within a 



projected circle with radius of ~1000AU. To carry the anal- 
ogy with the W3 region even further, we note that the ul- 
tracompact HII region/molecular h ot core W3(OH) lies at 
a similar distance from W3(H20) dWvrowski et al.l I1999I) as 
IRAS 05358+3543 mm2 and mm3 from mmla/b. The mea- 
sured expa nsion age of the W3(OH) UCHn region is only 
2300 years (Kawamura & Massor3ll998l) . It is well possible that 
IRAS 05358+3543 which now harbors a hypercompact Hn re- 
gion might evolve in to an UCHn region on a similar time scale, 
whereas mm2 and mm3 might evolve into regions resembling 
the present mmla/b [and W3(H20)] on similar time scales. 

The two sources mmla/b are located within one apparently 
common gas and dust core, and the projected separation of 
~1700 AU is w ithin the range o f observed distances for proto- 
binary sources dLaunhardtll2004l) . From the molecular line data 
presented in Leurini et al. (in prep.) we are not able to derive a 
clear velocity structure around the two cores which c ould be at- 
tributed to common rotation (e.g., Launhardt 2004). However, 
since we know about the complexity of the region and two 
molecular outflows emanating from the common mml core, it 
is not even a surprise that a coherent velocity signature is miss- 
ing. Each sub-source, mmla a nd mmlb, could drive one of two 
molecular outflows detected by B euther et al.l d200 2a). Based on 
the higher flux of mmla (Table [3]), its association with the com- 
pact cm continuum source, the Class n CH3OH maser and the 
main mid-infrared source, mmla is the main power-house of the 
region likely driving the more energetic large-scale collimated 
north-south outflow. 

Interestingly, of the four mm2 sub-sources, only mm2a is 
detected as a compact source in the highest resolution 1.2 mm 
image. This is indicative of mm2a being a genuine protostel- 
lar source whereas the emission features associated with mm2b, 
mm2c and mm2d could be more transient structures maybe 
caused by the underlying multiple outflow system. 

The source mm3 appears as the most simple object in the 
IRAS 05358+3543 massive star-forming region since it remains 
a single peak and shows a spectral index indicative of protostel- 
lar dust emission. The turnover of the SED at the shortest submm 
wavelengths indicates that mm3 is the coldest sub-source in this 
region. Hence, mm3 is a likely the youngest detected protostar 
in this system and may potentially drive one of the multiple out- 
flows identified in the region. 

While the CH3OH Class 11 maser can clearly be associated 
with mmla, none of the four H2O masers can be associated with 
any of the protostellar sources. Since H2O masers are believed 
to be produced in shock fronts they do not necessarily have to be 
located in close proximity to protostellar sources. 

4.2. Dust evolution and spectral indices 

While dust opacity indices /3 of 2 are usually found in the 
interstellar medium, lower values are possible in star-forming 
cores. For example, the dust opacit y per unit dust mass 
a:( 1.3m m) ~ 0.9cm 2 g _1 suggested by Ossenkorjf & Henning 
d 19941) for grains with thin ice mantles at gas densities of 
10 6 cirT 3 corresponds to /? ~ 1.4. Similarly, grain growth and 
dust coagulation in dense molecular cores decreases the values 
for j8 as well (e.g.. iBeckwifh et al.ll2000l lDraindl2006l) . A note 
of caution has to be added that the effect of a low spectral index 
a can also be caused by high optical depth from an underlying 
accretion disk and/or hypercompact Hn region, which can mimic 
low values of [5. Although it is difficult to differentiate between 
both scenarios, many studies favor the grain-growth scenario. 



The values /3 of 1.6 and 1.7 measured toward mml and 
mm3 from the lower-resolution dataset indicate that in these 
sub-sources the dust has probably not evolved much yet and the 
dust properties are still close to their properties in the interstellar 
medium. The SED of mm2 is too complex to derive a reasonable 
dust opacity index from the given dataset. 

The picture gets even more complicated going to the higher- 
spatial-resolution ratio map presented in Figure[6] Similar to the 
lower-resolution data, two of the protostellar objects (mml a and 
mm3) exhibit a rather normal spectral index of 3.6, correspond- 
ing to a dust opacity index /3 of 1.6. While the spectral index of 
2.7 for mmlb is already lower, it still is consistent with grain 
growth or high optical depth scenarios. 

Although the very low value of a ~ 2.3 toward the mm2a 
peak would still be consistent with grain growth in circumstellar 
disks or very high optical depths, the fact th at we see a similarly 
low spectral index on larger spatial scales (Willi ams et al.| [2004) 
indicates that additional physical processes may also come into 
play. Since we know of at least three molecular outflows ema- 
nating from the region, these outflows/jets are the most likely 
culprits to change the observed spectral indices. Thermal ion- 
ized wind models predict spectral indices between 2 and -0.1 
with typical sph erical winds exhibiting spectral indices a ~ 0.6 
(Reynolds 1986). Since these winds turn optically thin going to 
higher frequencies, their contribution to the SED is relatively 
larger at longer wavelength. Therefore, adding a wind compo- 
nent on top of the expected dust emission, a flattening of the 
spectral index is possible. Spectral indices below 2 as observed 
toward mm2c and its surrounding s can be due so lely to the out- 
flows within the region (see also Reynolds 1986). 

Even more puzzling are the large spectral indices > 4 
found in the vicinity of mm2b and north-west of mm2a. 
Similarly high spectral indices were observed toward the out- 
flo w region wi t hin th e LI 157 low-mass star-forming region 
by IGueth et al.1 (l2003h as well as toward some parts of the 
Orion molecular ridge and a cold molecular cloud near the 
Galactic Center (|Lis et al 1 1l998t iLis & Mentenl Il998h . While 
IGueth et ail (120031) attribute these high spectral indices to po- 
tentially smaller dust grains and/or chemical prope r ties (e.g., 
iKruegel & Siebenmorgenlll994l: IPollack et all fl994l ILis et all 
d!998l) and lLis & Mentenl (I 19981) find that their high-/? locations 
are at low temperatures. Therefore, they suggest that the high 
ft values can be caused by dust g rains covered with thick ice 
mantles (see also lAannestadll 975). Thi s picture is confirm ed by 
recent modeling of the dust opacity by iBoudet et al.l d2005l) in- 
dicating that p decreases with increasing temperature, and that (3 
can be as high as 2.5 (corresponding to a spectral index a = 4.5) 
at temperatures of 15 K. However, it is not obvious that the tem- 
peratures toward the high-alpha positions in Fig.|6]should be par- 
ticularly low, and thus it is questionable whether this picture can 
explain our observations. The fact that the positions of mm2b 
do not coincide between the 1 .2 mm and the 875 /im bands also 
indicates that dynamic processes like protostellar outflows can 
significantly alter the (sub)mm continuum emission. 

It is also interesting to compare the spectral indices we get 
at high spatial resolution with previous spectral index studies of 
the wh ole region with single -dish SCUBA observations at the 
JCMT Willia ms et all (|2004) derive a mean spectral index of 2 
increasing toward the (sub)mm peak to 2.6. This value is still 
significantly lower than the 3.6 we measure toward mml which 
is the position of the single-dish peak. Obviously, the single-dish 
data measure only some kind of weighted average of the values 
one resolves with interferometric observations. 



4.3. Density structure of massive star-forming cores 

The radial structure of the IRAS 05358+3543 dust continuum on 
large s patial scales p robed by bo lometer arrays h as been stud- 
ied by iBeuther eTail d2002bl) and lWilliams et alj d2005l) . These 
observations resulted in radial density profiles of the cluster- 
forming gas and dust clump with power-law distributions of 
p=l.9 and 1.5 for the 1.2 mm and 850 fj.m data, respectively. It 
is interesting that we find approximately the same power-law 
distributions in the uv-analysis of the individual sub-sources. 
Hence, the density structure of the large-scale cluster-forming 
clump appears to be similar to the small-scale structure of its 
fragmented star-forming sub-cores. Furthermore, the density 
distributions we find for the massive star-forming cores resem- 
ble th ose derived previously for low-mass star-forming regions 
(e.g.. iMotte & Andrei 1200 II) . The only other high-mass study 
we know of that analyzed the struct ure of massive sta r-forming 
cores with sub-arcsec resolution is bv lChen et alj d2006l) towards 
W3(H20) where they found density indices again close to 1.5. 
These results are consistent with the analyti c models and hydro- 
dyn amic simulation s of massive cores bv lMcKee&tanl {2003) 
and lKrumholz et al.l d2006l) . 

4.4. Limitations and future prospects 

Although these observations are among the highest-spatial- 
resolution multi-wavelength studies existing today, there are ob- 
viously still limitations to what we can resolve. The finest spatial 
resolution of ~ 0.5" we achieve at 1.2 mm corresponds at the 
given distance of 1800pc to a linear resolution of 900 AU, and 
we kn ow from young objects like W3IRS5 (e.g., Mege ath et al. 
2005) to more evolved regions like Orion (e.g.. iPreibisch et aL 
1999 ) that m u ltiplic ity at scales < 1000 AU can well be present. 
Kohl er et al.1 (|2006l) recently reported binary fractions for the 
Orion cluster of 23 + 10% for stars > 2M and 3.6 + 3.2% for 
stars < 2M Q . Even in IRAS 05358+3543 we resolve mml into a 
double-system (or triple system if including the MIR data) with 
spatial separations of ~1700AU ( 94.11 1. Therefore, all parame- 
ters we derive for the various (sub)mm sub-source have to be 
considered as potentially harboring multiple objects. This is par- 
ticularly important for the data at lower resolution (> 1"), where 
we know that they do not even resolve mmla and mmlb prop- 
erly. For example, as mentioned in 93.51 the radial density profile 
derived toward the peak position of mmla is actually the den- 
sity structure from the envelope encompassing mmla and mmlb. 
Similar arguments hold for the SED analysis and the mass and 
column density estimates. 

Nevertheless, it is likely that the expected multiple unre- 
solved sub-sources are gravitationally bound and that they form 
from the same gas and dust reservoir. Therefore, mass, col- 
umn density, temperature and density distribution of such en- 
tities are still the relevant parameters for the physical properties 
of the collapsing cores. This is consistent with recent 3D hy- 
drodynamically simulations showing that collapsing cores frag- 
ment less than previously expected, and that the most massive 
fragments are likely fo rming within the massive accretion disk 
(iRrumholz et al.ll2006l) which we do not resolve here anyway. 

The only way to observationally assess the fragmentation of 
young and cold massive star-forming cores in more detail is with 
higher spatial resolution and higher sensitivity. While the spatial 
resolution of the SMA and PdBI reaches today at best ~ 0.3" 
(less than a factor 2 better than the presented observations), the 
brightness sensitivity for these instruments at such scales is rela- 
tively poor. However, in only a few years from now, the Atacama 



Large Millimeter Array (ALMA) will be able to work regularly 
at 0.1" scales with good enough sensitivity to detect even lower- 
mass members of the evolving protoclusters. Then we will be 
able to derive the small scale-fragmentation processes in better 
detail than possible today. 



5. Summary 

We have observed the very young massive star-forming region 
IRAS 05358+3543 interferometrically with high spatial reso- 
lution in four frequency bands between 3.1mm and 438/im 
wavelength with the Plateau de Bure Interferometer and the 
Submillimeter Array. The observations reveal at least four sub- 
sources at the center of the region that are likely of protostellar 
nature. One of them, mmla, is associated with CH3OH Class 
11 maser emission, a cm wavelength hypercompact Hn region 
and a mid-infrared source. The center of two molecular out- 
flows - one of them the most collimated massive molecular out- 
flow observed so far in high-mass star formation - splits up 
into a double-source (mmla and mmlb) with a projected sep- 
aration of 1700 AU. With the given data, we cannot distinguish 
whether this projected double-source is really a proto-binary sys- 
tem or whether they are unbound. The additionally detected mid- 
infrared source ~1800 AU south-west of mmla indicates a rela- 
tively high protostellar density with at least three protostars in 
a projected circle with radius of ~1000AU. The main power 
house of the region is the sub-source mmla, but the other sub- 
sources are likely to power outflows as well. Furthermore, we de- 
tect additional (sub)mm peaks (mm2b, mm2c and mm2d) which 
may not necessarily be of protostellar nature but which could be 
caused by the underlying outflows and jets. H2O maser emission 
is found within the region of dust continuum emission but no 
feature is clearly associated with a protostellar (sub)mm contin- 
uum peak. This indicates that H2O maser emission is likely be 
produced in outflow shock fronts. 

This multi-frequency high-spatial-resolution (sub)mm con- 
tinuum dataset allows for the first time a spectral index analysis 
of individual sub-sources in a massive evolving cluster that con- 
tain either individual protostars or unresolved, likely gravitation- 
ally bound multiple systems. The spectral energy distribution of 
the main power house mml resembles a typical SED one ex- 
pects for a massive star-forming core: a hypercompact Hn region 
from a 13 M Q Bl star is detected at cm wavelength whereas the 
mm to submm spectrum is dominated by a cold dust component 
with a spectral index of 3.6. The 438/im data are of particular 
interest because they fit well into the dust spectrum derived pre- 
viously from the longer wavelength data. The SED of the sub- 
source mm3 shows a different behavior at this short wavelength: 
although we can only derive upper limits at 438 //m for mm3, 
this upper limit is significantly below the power-law spectrum 
that one can fit from the lower-frequency data and that would be 
expected for a normal Rayleigh-Jeans-like part of the SED. Such 
a lower flux measurement can be explained by even colder dust 
emission where the peak of the SED shifts to longer wavelength 
and hence the 438 fim point is already closer to this peak and 
thus not on the Rayleigh-Jeans tail of the Planck-function any- 
more. Trying various grey-body functions, we estimate an upper 
limit for the dust temperature of 20 K. 

A radial analysis of the data in the uv-domain reveals that 
the density power-law indices p vary for the different sub- 
sources between 1 .5 and 2, similar to observed density distribu- 
tions of low-mass star-forming cores. Furthermore, these power- 
law indices resemble those previously found for the larger-scale 



cluster-forming clump, indicative of similar density distributions 
on the high-mass star- and cluster-forming scale. 

Although the discussed sub-sources may still harbor mul- 
tiple objects, the derived physical parameters should repre- 
sent the properties of individual collapsing massive cores well. 
Fragmentation on even smaller spatial scales (e.g., within the 
massive accretion disks) will likely only be observationally ac- 
cessible with the advent of ALMA in a few years from now. 
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Fig. 4. (Sub)mm continuum images in the four observed bands 
smoothed to the spatial resolution of 1.9" x 1.4" of the low- 
est resolution image at 3.1mm. The contour levels are from 
10 to 90% (step 10%) of the peak emissions with peak fluxes 
of 10.3, 91.8, 270.8, and 3569mJy beam 1 for the 3.1mm, 
1.2 mm, 875, and 438 /im bands, respectively. The 438 ^m im- 
age shown here is not primary-beam corrected yet, the size 
of the primary beam is indicated by the large half-circle in 
the bottom panel. The left and right vertical lines mark the 
positions of mml and mm3. In contrast, the peak position 
of mm2 shifts with wavelengths in these lower-resolution im- 
ages which is indicated by the three middle vertical lines. 
However, the full circles mark the mm2a peak position in the 
corresponding high-resolution images (Fig. |2). The files for 
this figure are also available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat? J/A+A/. 
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Fig. 6. Ratio map of the 875 //m versus the 1.2 mm data. Both 
datasets were smoothed to a common spatial resolution of 1 .2" x 
0.8". To better align the two maps before producing the ratio, the 
875 //m was shifted by 0.18" to the west. The top panel shows 
the two maps at the same spatial resolution. The bottom panel 
then presents the ratio map in color scale with the 875 yum data at 
original resolution overlayed in contours. The synthesized beam 
of the data used for the ratio map is shown at the bottom-right 
of the top panel, and the bottom-right of the bottom panel shows 
the synthesized beam of the 875 pm map at original resolution. 
The files for this figure are also available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79. 128. 5) 
or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat? J/A+A/. 
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Fig. 7. Plots of amplitude versus uv -distance for the three main 
sub-sources in IRAS 05358+3543 at three observed frequencies. 
Results from power-law fits are shown as well. Although the top 
panel was produced centered on mm la, the data <250k/l are 
from scales combining mmla and mmlb. Therefore it resembles 
the intensity and density structure of this multiple system. 



